Kondo memory in driven strongly-correlated quantum dots 
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We investigate the real-time current response of strongly-correlated quantum dot systems under 
sinusoidal driving voltages. By means of an accurate hierarchical equations of motion approach, 
we demonstrate the presence of prominent memory effects induced by the Kondo resonance on the 
real-time current response. These memory effects appear as distinctive hysteresis line shapes and 
self-crossing features in the dynamic current-voltage characteristics, with concomitant excitation of 
odd-number overtones. They emerge as a cooperative effect of quantum coherence - due to inductive 
behavior - and electron correlations - due to the Kondo resonance. We also show the suppression 
of memory effects and the transition to classical behavior as a function of temperature. All these 
phenomena can be observed in experiments and may lead to novel quantum memory applications. 

PACS numbers: 72.15.Qm, 73.63.-b, 85.35.Be 



CD 



I 

o 
o 



> 

o 
o 

(N 

o 
m 



X 



Memory of physical systems is quite a common phe- 
nomenon [l| . It simply means that the state of a system 
at a given time - driven by an external input - depends 
strongly on its history, at least within a certain range of 
parameters of the external drive, such as its frequency 
or amplitude [1, 2]. These memory features give rise to 
unconventional properties and novel functionalities that 
can be used in actual device applications. 

Although there has been a recent crescendo of inter- 
est in memory effects in resistive, capacitive and induc- 
tive systems 3], less work has been devoted to memory 
effects in strongly-correlated quantum systems. These 
include the single- molecule magnets [J, |5| , metal-oxide- 
semiconductor thin films [GJ, |7|, nanoparticle assemblies 
Isj, and bulk Kondo insulators [9J. In this respect, quan- 
tum dots (QDs) coupled to electron reservoirs are ideal 
platforms for experimental and theoretical studies rele- 
vant to quantum computation and quantum information 
lOl lll |. This is because the electron or spin state of a 
QD can be manipulated conveniently by external fields, 
with the quantum coherence largely preserved. In partic- 
ular, at low temperatures, formation of the Kondo singlet 
state 12h15| opens up additional channels for electron 
conduction. However, it is not at all obvious if new fea- 
tures would emerge due to the interplay between quan- 
tum coherence and electron correlations under driving 
conditions, and if the Kondo phenomena would influence 
significantly the memory of the QD system. 

In this letter we show that Kondo resonances indeed in- 
duce prominent memory effects when a QD is driven by 
an external periodic voltage. These features appear as 
well-defined hysteresis line shapes and self-crossing fea- 
tures in the dynamic current-voltage characteristics with 
excitation of odd-number overtones. They are a direct 
consequence of quantum coherence and electron correla- 
tions and are suppressed at relatively high temperatures 



where a classical behavior is recovered. All these ef- 
fects can be probed at specific frequencies that are within 
reach of experimental verification. 

The QD system of interest is illustrated with a single- 
impurity Anderson model [l6|, i?dot = ed("-t + "-i) + 



Ufi^fi^. Here, 



and a\ (Ss) creates (anni- 



hilates) a spin-s electron on the dot level of energy e^, 
while U is the on-dot electron-electron (e-e) interaction 
strength. The total Hamiltonian is 
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J2ak takdldak + ii-C., for the nonintcracting leads and the 
dot- lead couplings, respectively. Here, dj^i, (dak) is the 
creation (annihilation) operator for the a-lead state |fc) of 
energy eak , and tak is the coupling strength between the 
dot level and |fc). For numerical convenience, identical 
leads are considered, and their hybridization functions 
assume a Lorentzian form, Aq,(cj) = tt^^ tafci*^. (5(w — 
eak) = AW'^/2[{uj ~ ^XaY + W^], where A is the effective 
dot-lead coupling strength, W is the band width, and /Xq, 
is the chemical potential of the a-lead. We prepare the 
initial total system at equilibrium, where ^a = M°'^ = 0- 
An ac voltage is applied to the left (L) and right (R) leads 
from the time to = 0, i.e., Vi,{t) = — VR(t) — Vq sin{VLt), 
with Vb and ^ being the voltage amplitude and frequency, 
respectively. Upon switching on the voltage, the QD is 
driven out of equilibrium, and the time-dependent cur- 
rent flowing into the a-lead Ia{t) is computed. 

To proceed, a choice of method for calculating current 
must be made that accurately reproduces the nonequilib- 
rium electronic dynamics of this strongly correlated elec- 
tron system. A number of theoretical approaches have 
been developed to tackle this problem, such as the time- 
dependent numerical renormalization group (NRG) ap- 
proach [17l | , the time-dependent density matrix renormal- 



ization group approach 18|, 111 
matic quantum Monte Carlo 2 



the real-time diagram- 



21| , the KadanofF-Baym 



time-dependent Green's function approach 22|, |23| , the 
time-dependent density-functional theory [24 [ , and quan- 
tum master equation approaches [5|, I25l427l | . For instance, 
the time scale for the buildup of Kondo resonance in 
a strongly correlated QD has been studied by using a 
Green's function technique, where the e-e interaction 
is accounted for within the noncrossing approximation 

We choose the hierarchical equations of motion 
(HEOM) approach developed in recent years by Yan 



M^- 



and coworkers |3Cll - l32| . This is a unified, real-time ap- 
proach for a wide range of equilibrium and nonequilib- 
rium, static and dynamic properties of a general open 
quantum system |32h35| . The basic variables in HEOM 
are the reduced density matrix of the open system and a 
set of auxiliary density matrices. The HEOM approach 
has been employed to study, for instance, the dynamic 
Coulomb blockade and dynamic Kondo phenomena in 
QDs |36|, |37|. To close the coupled equations, the hi- 
erarchy needs to be truncated at a certain level L. In 
principle, the exact solution is guaranteed for L -> oo. 
In practice, the results usually converge rapidly with in- 
creasing L at finite temperatures. Once the convergence 
is achieved, the numerical outcome is considered to be 
quantitatively accurate. For instance, it has been demon- 
strated that the HEOM approach achieves the same level 
of accuracy as the latest NRG method for the prediction 
of various dynamical properties at equilibrium [32| , and 
it has reproduced [33[ the exact numerical solution for 
the time-dependent current response of a noninteract- 
ing QD to a step-function bias voltage 38|,|39|. Here, we 
solve the HEOM dynamics described by Hamiltonian ([1]) , 
at the L — 4 truncation level of hierarchy 36l ISTj (see 
also Supplemental Material [40[). We have checked that 
all results presented in this letter converge quantitatively 
with respect to this truncation. 

Figure [1] depicts the dynamic I-V characteristics of a 
QD which possesses the electron-hole symmetry {U = 
—2ed) subject to an ac voltage of various frequencies il. 
To compare with the energetic parameters of the QD, 51 is 
represented by the harmonic energy associated with the 
driving voltage. In particular, J7 = 0.01 meV corresponds 
to a period of 0.4 ns. Although high, these GHz time- 
domain voltage manipulations and current measurements 
have been realized experimentally 4l|-|43|. The calcula- 
tion results indicate that the electron-hole symmetry is 
preserved under the ac voltage, and the dot level is al- 
ways half-filled. Consequently, the displacement current 
is zero, and the Kirchhoff 's current law holds at any time, 
i.e., /(t) = /R(i) = -/L(i) i. 

As shown in Fig.[lja), at a low il, the dynamic I-V 
characteristics is close to the steady-state curve, since 
the electrons have sufficient time to redistribute to catch 
up with the variation of voltage. In contrast, at a high 
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FIG. 1: (Color online). Current versus voltage for a QD 
driven by an ac voltage of amplitude Vo — 0.1 mV and fre- 
quency SI. The arrows indicate the circulating directions of 
the loops. The green dashed curve in (a) depicts the steady- 
state I-V characteristics. The insets plot the corresponding 
current (in nA) versus time (in ps). The other parameters (in 
meV) are U = -2ed = 0.4, A = 0.1, W = 2, and T = 0.02. 



n, the dynamic I-V curve forms an ellipse; see Fig.[TJd). 
The ellipse shape originates from a phase difference be- 
tween the current and voltage. In the case of Fig.[Ud), 
the ellipse is traversed in the counterclockwise direction 
(voltage leads current), indicating a prominent inductive 
behavior. It has been demonstrated that in the linear 
response regime a two-terminal device can be mapped 
onto a classical circuit, where a resistor-capacitor branch 
and a resistor-inductor branch are connected in paral- 
4414461. Here, the inductance is associated with the 



lei 

time for an electron to dwell on the device [4^, |47[ and 
is due to the fact that the dot energy is lower than the 
equilibrium energy, as required in order to have a Kondo 
resonance. At an intermediate fl, the dynamic I-V curve 
forms a hysteresis loop with self-crossing in the first and 
third quadrants; see Fig.fljb). The hysteresis behavior 
highlights the nontrivial memory effects on the real-time 
electron dynamics. To understand further these memory 
effects, we vary both the coupling strength U and the 
temperature T. 

Figure [5] plots the dynamic I-V curves for a series 
of symmetric QDs of different U under the same driv- 
ing voltage and temperature T. The hysteresis behav- 
ior of Fig.[ljb) is constantly observed. It is interest- 
ing to see that the I-V curve shape varies in a non- 
monotonic manner with respect to U, and the hysteresis 
behavior is most outstanding at around U = 1.2meV. 
To understand this trend, we calculate the steady-state 
differential conductances {dl/dV) of the QD under the 
bias voltage Vl = — Vr, = V. The resulting dl/dV 
are then scaled by their respective values at T^ = 0. 
These curves are shown in Fig.[21[d). The differential 
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FIG. 2: (Color online). Current versus voltage for the QDs 
of (a) U = 0.4, (b) U = 1.2, and (c) U = 2meV, under the 
antisymmetric ac voltage of Vq = 0.05 mV and Q, = 0.03 meV. 
Panel (d) shows the steady-state dl /dV versus the bias volt- 
age V = Vl = —Vn, where the displayed data are scaled by 
the value of dl/dV at V^ = 0. The other parameters (in meV) 
are ea = -U/2, A = 0.1, M/ = 2, and T = 0.01. 



conductance clearly decays nonlinearly with increasing 
voltage, and U = 1.2 meV corresponds to the steepest 
descent of dl/dV. The maximal conductance at t/ = 
is due to the formation of Kondo states. Under a bias 
voltage, the Kondo channel is blocked by the gap be- 
tween the chemical potentials in the two leads, leading 
to the drastic decay in conductance with increasing V. 
This is confirmed by the evolution of the Kondo spec- 
tral signature versus voltage (see Supplemental Mate- 
rial [40|)- For the QD studied in Fig. [21 the Kondo tem- 
peratures Tk = (f/A/2)i/2exp[-7rC//8A + 7rA/2t/] are 
0.044, 2.5 X 10~3 and 1.3 x lO^^meV for U ^ 0.4, 1.2 
and 2 meV, respectively. Therefore, under a fixed tem- 
perature, the Kondo resonance would lead to a steeper 
decay in dl/dV with a smaller U. On the other hand, 
the dot level of energy e^ = ~U/2 becomes closer to 
fi°'^ as U diminishes. This provides another conduction 
channel with width A, and its associated conductance 
decays much more slowly than the Kondo channel since 
Tfj- ^ A. With the presence of both conduction chan- 
nels, the steepest decay in dl /dV occurs at an interme- 
diate value of U. 

The shape of the dynamic I-V characteristics can now 
be interpreted as follows. In the limit of adiabatic driving 
(ri — >■ 0), the evolution of the current response follows ex- 
actly the steady-state dl/dV-V curve. The reduction in 
dl/dV versus V leads to a concave I-V curve at F > 0; 
see the green dotted line in Fig.[lja). At a finite J7, the 
inductive feature starts to emerge, which forces the cur- 
rent to lag behind the driving voltage. For instance, at 
an instant when the voltage turns from negative to pos- 
itive, the current has not yet reached its turning point 



- it remains negative but accelerates toward the volt- 
age. Consequently, the dynamic I-V curve has a convex 
curvature at F = and / < 0. The convexity may 
be inverted into concavity provided that the increase in 
voltage captures the drastic decay in conductance. This 
requires a relatively low fi, with which the variation of 
"transient" conductance does not deviate much from the 
steady-state dl/dV-V curve shown in Fig.[2][d). Such a 
convex-concave curve segment, together with its mirror- 
image counterpart upon voltage decrease, results in a 
self-crossing hysteresis loop as depicted in Fig.[ljb) and 
Fig.[2][a)-(c). At a high ft, there is not enough time for 
the buildup of Kondo states, and thus the dl/dV-V curve 
is effectively smoothed. In this case, the inductive fea- 
ture becomes a dominant factor, and the dynamic I-V 
curve displays a simple ellipse shape; see Fig.[T][d). 

Based on the above analysis, we can conclude that the 
hysteresis behavior of the ac I-V characteristics origi- 
nates from a cooperative interplay between the quantum 
coherence (the inductive feature) and strong electron cor- 
relation (the Kondo resonance). Since both quantum 
coherence and correlation are significantly infiuenced by 
thermal fluctuations, the associated memory effects are 
expected to depend sensitively on the temperature. 

In Fig.[31[a) we show the evolution of the dynamic I- 
V curve with increasing temperature. Here, in order to 
accentuate the nonlinear current response, a somewhat 
larger voltage amplitude {Vq = 0.1 meV) is adopted. At 
low temperature (T = 0.01 me V), the ac I-V curve dis- 
plays a remarkable hysteresis behavior, with the emer- 
gence of multiple self-crossing points. Such a convo- 
luted line shape is a strong indication of complex mem- 
ory effects. With the temperature rising by as little as 
0.01 meV, the hysteresis behavior of Fig.[2l^b) is recov- 
ered, with only one crossing in the first or third quadrant. 
This apparent change in the hysteresis pattern should be 
ascribed to the suppression of the Kondo resonance. In- 
deed, as shown clearly in Fig.jSjc), the zero-bias conduc- 
tance decays drastically with a minor increase in T from 
0.01 meV. Meanwhile, increasing T leads to the progres- 
sive narrowing of the I-V hysteresis loop, and the loop 
almost merges into a single line at T = 0.1 meV. This 
means that the inductive feature becomes irrelevant at 
a high temperature, and the QD behaves like a classical 
resistor. This is because the conducting electrons have a 
wide energy distribution due to thermal excitations, and 
the phase difference between the current and voltage is 
largely randomized. Due to the same reason, the steady- 
state dl/dV varies rather smoothly as a function of V at 
high T, leading to an overall Ohmic conductance. 

The temperature-sensitive memory is further explored 
through a frequency analysis of the response current. The 
current spectrum for a single-lead QD under sinusoidal 
voltages has been studied in Ref.|45|, where the response 
current signals concentrate on the exact multiples of driv- 
ing frequency 51. For the two-terminal QD investigated 
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FIG. 3: (Color online), (a) Dynamic I-V characteristics, (b) 
real-time current response, (c) steady-state dl /dV versus V, 
and (d) imaginary part of current spectrum of a QD with 
U — —2ed = 1.2 meV at various temperatures. The inset of 
(a) magnifies the self-crossing of ac I-V curve near V = — Vo 
at r = 0.01 meV, and that of (b) displays the full peaks of 
Im[7(a;)] at w = il. The other parameters (in meV) are A = 
0.1, W = 2, l/o = 0.1, and Q. = 0.03. 



here, the variation of current spectrum versus temper- 
ature is depicted in Fig.[31Jd). It is found that the re- 
sponse current arises predominantly at w = (2n -I- 1)0 
with n G Z, while the even overtones remain unexcited. 
This is because the ac voltages have the left-right anti- 
symmetry, so that the response current would flow in 
the reverse direction if the half-period were taken as 
the initial time (ip — 7r/f2). Therefore, the relation 
liu}) = —I{w)e'^'^'"/^ must be satisfied. Overtones of 
all symmetries can be found by breakin g th e symmetric 
coupling between the QD and the leads [48]. As T rises 
from 0.01 to 0.1 meV, the change in current amplitude at 
the fundamental frequency is rather minor. In contrast, 
the overtone responses vanish almost completely. This 
highlights the dynamic Kondo-assisted electron conduc- 
tion, which requires a relatively low excitation energy. 
At T < 0.01 meV, the Kondo memory effects are ex- 
pected to be even more prominent. However, to achieve 
convergent HEOM results at a lower T, a higher trun- 
cation level L is required, which is numerically too de- 
manding with the present computational resources at our 
disposal. At T > 0.1 meV, the Kondo resonance is ab- 
sent, and the dot level is away from the resonance region 
(|ed| > Vq). Consequently, the system falls in the linear 
response regime, and the corresponding real-time current 
becomes synchronized with the voltage; see Fig.[3l[b). 

We finally present some considerations on the exper- 
imental verification of our predictions. In realistic ex- 
perimental setups, a QD may involve more than one en- 
ergy level and more complicated forms of e-e interaction. 
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FIG. 4: (Color online). Dynamic I-V characteristics of QDs 
with different e^. The displayed current is I{t) = [Init) — 
Ih(t)]/2. The other parameters (in meV) are A = 0.1, U = 
1.2, T = 0.01, W ^2, Vo = 0.05, and Q, = 0.03. 



Therefore, the memory effects will influence the real-time 
dynamics in a rather complex way. However, since the 
Kondo resonance is tightly tied to the lead chemical po- 
tential, the Kondo memory may be identified distinctly 
from other memory effects. For instance, consider the 
scenario where the dot-level e^ is shifted by a gate volt- 
age, and the ac I-V curves at different e^ are shown in 
Fig.m At U ^ — Scd, the electron- hole symmetry is bro- 
ken, and transient charging takes place at the QD. To 
concentrate on the net conduction current, the displace- 
ment component is taken care of by symmetrizing the 
left- and right-lead currents J5J. As clearly indicated in 
Fig-Hl the Kondo-resonance-induced hysteresis behavior 
is continually preserved with e^ shifted by as much as 
0.4 meV. Since any single-electron resonance associated 
with a dot energy level would have been switched off 
under such a large gate voltage. Fig. 2] accentuates the 
robustness of Kondo memory effects. 

To conclude, we have demonstrated the presence of 
Kondo-resonance-induced memory effects in the real- 
time electronic dynamics of strongly correlated QDs, 
which may be observed in experiments. The predicted 
hysteresis line shape and self-crossing features of the ac 
I-V characteristics highlight the remarkable interplay be- 
tween quantum coherence and quantum electron correla- 
tions. These effects may lead to novel applications of QDs 
in the Kondo regime, ranging from machine learning [49| 
to massively-parallel computation and information pro- 
cessing [50,]. 
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